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The global minima of clusters bound by a Dzugutov potential form non-compact polytetrahedral 
clusters mainly composed of interpenetrating and face-sharing 13-atom icosahedra. As the size 
increases, these icosahedral units first form linear arrays, then two-dimensional rings, then three- 
dimensional networks. Characterization of the energy landscapes of these clusters shows that 
they are particularly rough and generally exhibit a multiple-funnel topography. These results 
provide new insights into the structure and dynamics of bulk supercooled Dzugutov liquids and 
the form of the bulk phase diagram. 
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I. INTRODUCTION 

One reason for the continued growth in cluster sci- 
ence is the link that these finite systems provide be- 
tween atoms and molecules and bulk systems. Indeed, 
much cluster research has focussed on the size evolution 
of cluster properties in order to understand how bulk be- 
haviour, such as phase transitions and crystalline struc- 
ture, emerge with increasing size. Furthermore, from this 
research new insights into bulk behaviour can often be 
obtained. 

For example, the structure of clusters has often been 
invoked when considering the local structure of bulk liq- 
uids and glasses. Frank may have been the first to make 
this connection^ He noted that the 13-atom icosahedron 
has a lower energy than a close-packed cluster for the 
Lennard-Jones (LJ) potential. From this result he in- 
ferred that the most favoured structure of liquids would 
involve icosahedral local ordering, thus making homoge- 
neous nucleation of a crystal from a liquid difficult, be- 
cause crystal formation would require a change in the 
local order a Since then convincing evidence that simple, 
liquids have polytetrahedraH order has been obtained,cl 
and so the use of clusters as models of local order in liq- 
uids has become commonfl Furthermore, clusters with 
five-fold symmetry are ofteujised as building blocks for 
quasicrystalline structuresffijtra and metallic crystals.!^ 

Much of the computational work on supercooled liq- 
uids and glasses has focussed on the behaviour exhibited 
by systems interacting according to simple model poten- 
tials. One of the most commonly studied systems em- 
ploys a binary LJ potential,0 where the binary nature of 
the system hinders the crystallization that readily occurs 
for a single component LJ liquid. However, recent evi- 
dence suggests that the low energy configurations for this, 
system involve some demixing and crystalline ordering.E3 



Furthermore, although this system is a good generic glass 
former, it does not capture all of the features associated 
with metallic glasses, for which particularly strong local 
polytetrahedral and icosahedral ordering occurs. 

An alternative approach is to use a single compo- 
nent system with a potential that is designed to prevent 
crystallization. This line of enquiry was pioneered by 
Dzugutov, who introduced a potential with a local maxi- 
mum at approximately \J~2 times the equilibrium pair^dis-, 
tance in order to disfavour close-packed structures. 
This maximum also bears a resemblance to the first of 
the Friedel oscillations that can occur for metal poten- 
tials. Systems interacting according to this Dzugutov 
potential have been found to be good glass-formers, and 
to exhibit a first sharp diffraction-eeak and a split sec- 
ond peak in the structure factorjij which are common 
features of metallic glasses. One other interesting fea- 
ture of this potential is that under certain conditions it 
can produce a dodecagonal quasicrystal.ES 

Our aim here is to further understand the structural 
implications of the Dzugutov potential by identifying the 
characteristic structural motifs associated with the global 
minima of Dzugutov clusters. Furthermore, we exam- 
ine the consequences of the potential for the topogra- 
phy of the energy landscape [or potential energy surface 
(PES)] of these clusters. As well as the insights into 
bulk behaviour that these results provide, they will also 
be potentially helpful in understanding the structures of 
clusters that favour polytetrahedral structure. Recent 
results on smalL-cobalt clusters indicate that they are 
polytetrahedraljLZI although the detailed structure is un- 
known. Furthermore, one would expect clusters of alloys 
that adopt a Frank-Kasper structureMESIij in bulk to 
show polytetrahedral order in clusters. 
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FIG. 1: A comparison of the Dzugutov (Dz) and 
Lennard-Jones (LJ) potentials. The unit of energy is the 
pair well depth. The vertical lines show the positions of 
the pair distances in the 13-atom icosahedron and their 
heights are proportional to the numbers of pairs with 
that distance. 



II. METHODS 



A Potentials 



The Dzugutov pair potential is^° 1 
V(r) = A{r- rn - B) exp 
B exp 



9(a-r) + 



r — b 



6(6 -r), 



(1) 



where Q(x), the Heaviside step function, is for x < 
and 1 otherwise, and the parameters in the potential have 
the values 



,4 = 5.82 a = 1.87 

5 = 1.28 6=1.94 

c=1.10 d = 0.27 

m = 16. 



(2) 



The total potential energy of a cluster of Dzugutov atoms 



tween atoms i and j 
mum at r, 



is then E = J2i<i ^( r y)> where is the distance be- 
The pair potential has a mini- 
q— 1.130 of depth e=0.581 and a maximum at 



r max =l. 628=1.44 r oq of height 0.460 (0.791e). 

In Figure [l] we compare the Dzugutov and the LJ po- 
tentials. The curvature at the bottom of the Dzugutov 
well in the above units is virtually the same as that of the 
LJ potential in reduced units. However, when the unit 
of energy for both potentials is the pair well depth, as in 
Figure |l|, the curvature at the minimum of the Dzugutov 
potential is somewhat larger, i.e. the Dzugutov potential 
well is narrower. 



This effect can be quantified. In units of pair well 
depth and equilibrium pair separation, the curvature at 
the minimum of the Morse potential 
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is 2p and of the generalized LJ potential 
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V(r) = 4e 
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is 2n 2 . These results allow one to calculate an effective n 
(or p) value. For the Dzugutov potential n^|=pj^|=7.52, 
as compared to n LJ =6. This comparison will be impor- 
tant to our interpretation of the structures p£|thcse clus- 
ters since, as illustrated for Morse clusters, EHEil the well 
width is a key factor in determining a cluster's structure. 

For a paiiHpotential we can partition the energy into 
three terms:c2l 



E = -n nn e + E. 



strain 



E n 



(5) 



where e is the pair well depth, n nn is the number of near- 
est neighbours (two atoms are defined as nearest neigh- 
bours if r-ij < ro), 



^strain = ^ ( V ( r v) + € ) aild ( 6 ) 

(7) 



i<3,rij<ro 



£„„„= v ^)- 

i<j,rij>r 



We choose ro=1.378 where V(ro)=0. The first term in 
Equation (^|) is the ideal pair energy if all n nn nearest- 
neighbour pairs lie exactly at the equilibrium pair dis- 
tance. The strain energy, -Entrain, is the energetic penalty 
for the deviation of nearest-neighbour distances from the 
equilibrium pair distance, and -Ennn is the contribution 
to the energy from non-nearest neighbours. For the 
Dzugutov potential E nnn corresponds to the energetic 
penalty for distances lying near to y/2 r cq (more precisely 
in the range r*o < r < b). The global minimum of a clus- 
ter represents the best balance between maximizing n nn 
whilst minimizing both the strain energy and £nnn- 

The big difference between the Dzugutov potential and 
most other pair potentials is the local maximum at r max 
and the consequent positive value of the term E nnn . This 
term disfavours distances near to ^/2r cq that are as- 
sociated with the octahedral interstices in close-packed 
structures—For clusters this term also disfavours Marks 
decahedraca and MackayE3 icosahedra with more than 
one shell, because the strained face-centred cubic (fee) 
tetrahedra that make up these structures contain octa- 
hedral interstices. However, as Figure [l] shows, the dis- 
tances in the 13-atom icosahedron are able to avoid this 
maximum. Therefore, E nnn favours polytetrahedral clus- 
ters. For larger polytetrahedral clusters the value of E nnn 
is somewhat correlated with -Entrain- Strain results from 



3 



a dispersion of nearest-neighbour distances, and so is also 
likely to lead to a greater dispersion in the second nearest- 
neighbour distances and potentially more distances that 
stray towards the maximum at r max . 

We have previously analysed the roles of n nn and 
-Entrain fop, .other pair potentials, such as the Morse 
potential. E9EII A large value of n nn is achieved by hav- 
ing a compact spherical shape with a large proportion of 
surfaces with high coordination numbers, e.g. fee {111} 
faces, ^strain disfavours strained structures and the mag- 
nitude of this term increases as the potential well be- 
comes narrower. The strain for close-packed structures 
can be zero, but increases along the series: decahedral, 
icosahedral and polytetrahedral. Thus, the strain associ- 
ated with compact polytetrahedral structures can be only 
accommodated at low n e g values. Therefore, given the 
value of n®g , compact polytetrahedral structures seem 
unlikely and it is not self-evident what types of structure 
are best able to satisfy these competing demands. 



tential exhibit highly non-spherical geometries it was nec- 
essary to enlarger-the container to twice the radius used 
in previous work.u Atoms with pair energies of zero were 
therefore shifted towards the origin through a distance of 
one unit on every call of the potential subroutine. 

C Searching the potential energy surface 

To map the PES topography of some of the Dzugutov 
clusters we used the s ame . m ethods as these we have pre- 
viously applied to LJl23'E§Ea and MorseS3 clusters. We 
thereby obtain large samples of connected minima and 
transition states that provide good representations of the 
low-energy regions of the PES. The approach-anvolves 
repeated applications of eigenvector-following to hnd 
new transition states and thp-.minima they connect, as 
described in detail elsewhere£3 



B Global Optimization 

To obtain putative global minima-, for the Dzugutov 
clusters we used tie basin-hoppingc-3 (or Monte Carlo 
plus minimizationEj) algorithm, which has proved very 
useful in locating cluster global minima.E3 In this algo- 
rithm a constant temperature Metropolis Monte Carlo 
simulation is performed with the acceptance criterion 
based not upon the energy of the new configuration, but 
upon the energy of the minimum obtained by a local 
minimization from that configuration. Details of the ba- 
sic approach employed in the present work may be found 
elsewhere, l-3 

Dzugutov clusters are particularly challenging cases for 
global optimization, and so for each cluster size five runs 
of 100000 steps each starting from a random configu- 
ration were performed. Short runs of 10 000 steps were 
also performed for each cluster starting from the lowest- 
energy minima found for the size above and below with 
an atom removed and added, respectively. In the lat- 
ter runs the seed geometries were fixed for the first 2 000 
steps. These short runs were particularly important at 
larger sizes, and were repeated until there were no fur- 
ther improvements for any size. In each run a constant 
temperature of 0.4 in reduced units was employed, with 
a fixed maximum step size for Cartesian displacements 
of 0.36. Random angular moves were performed for the 
most weakly bound atom if its pair energy exceeded 3/5 
of the lowest pair energy for the structure in question. 
(The pair energy of atom a is defined as J2i=£ a V( r ia)-) 

Problems arise in the energy minimizations if atoms 
evaporate from the cluster: there are no attractive forces 
on an atom lacking any neighbours within a distance of 
r max . As in previous work we constrained the atoms in 
a sphere to limit the maximum interatomic distances. 
However, since some global minima of the Dzugutov po- 



III. RESULTS 



A Global Minima 



The lowest-energy minima that we found are given in 
Table | up to A=100; these results will alscLbe made 
available in the Cambridge Cluster Database.!!^ For most 
clusters up to N ~ 60 these structures were obtained 
from basin-hopping runs that were started from a random 
configuration. However, above this size we often found 
that the lowest-energy minimum was only located in the 
runs seeded with a configuration derived from the best 
structure at a smaller or a larger size. Therefore, it would 
not be surprising if lower-energy minima exist for some 
of the larger clusters. However, any such minima would 
not be expected to change the observed structural motifs, 
which are the main focus of the present section. 

The energies of the putative global minima are plotted 
in Figure relative to the smooth function, i? avo , which 
represents the average energy of the global minima, so 
that particularly stable sizes are readily apparent. A rep- 
resentative set of global minima is depicted in Figure |^; 
most of these structures correspond either to a minimum 
in Figure |^ or a high symmetry configuration. 

As noted earlier, one consequence of the maximum in 
the potential is that polytetrahedral structures are likely 
to be favoured, and for small sizes the clusters adopt 
the expected structures. The 13-atom global minimum 
is the icosahedron (Figure ||) and the global minima at 
smaller sizes lie on a polytetrahedral growth sequence 
to this structure. This sequence is similar to the global 
minima for the L J potential^ except that the Dzg global 
minimum is the bicapped tetrahedron, rather than the 
octahedron. 

There are two possible overlayers for growth on the 
13-atom icosahedron that are frequently observed. The 
first involves adding atoms to the faces and vertices of 
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the icosahedron and is termed the anti-Mackay overlayer. 
This overlayer continues the polytetrahedral packing of 
the icosahedron and leads to a series of interpenetrat- 
ing icosahedra. The overlayer is completed at a 45-atom 
rhombic tricontahedron, which is an icosahedron of inter- 
penetrating icosahedra. The second possible overlayer in- 
volves adding atoms to the edges and vertices of the icosa- 
hedron, and is termed the Mackay overlayer because it 
leads to the 55-atom Mackay icosahedron. This overlayer 
introduces octahedral interstices and so is disfavoured by 
the Dzugutov potential. No Dz n global minima have this 
structure. 

The global minima above TV=13 initially have an anti- 
Mackay structure, for example, the double icosahedron 
at iV=19 (Figure |). This trend continues up to TV=22. 
For TV=23 the usual anti-Mackay structure is a trian- 
gular arrangement of three interpenetrating icosahedra. 
However the DZ23 global minimum consists of two face- 
sharing icosahedra (Figure ||), because the anti-Mackay 
structure has a larger strain energy (4.810 as opposed 
to 3.138 for the global minimum) and a slightly higher 



£ nnn , making it 0.156 higher in energy, even though it 
has three more nearest-neighbour contacts. 

This trend is continued for larger sizes. The strain 
associated with compact polytetrahedral structures can- 
not be accommodated without too great an energetic 
penalty because of the value of njpff . Instead, polytetrahe- 
dral structures with extended configurations are favoured 
even though they have fewer nearest neighbours. These 
structures involve a combination of face-sharing and in- 
terpenetrating icosahedra (Figure ||). At TV=25 a bent 
arrangement of three interpenetrating icosahedra is the 
global minimum (it is just lower in energy than the lin- 
ear triple icosahedron). At TV=31 and 37 the global min- 
ima are linear quadruple and quintuple icosahedra. At 
TV=33 the global minimum is a linear arrangement of 
three face-sharing icosahedra. At TV=29 the global mini- 
mum involves an icosahedron sharing a face with a double 
icosahedron and at TV=35 it is two face-sharing double 
icosahedra. 

Above N=37 the global minima become two- 
dimensional in character. The 38-atom global minimum 
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FIG. 2: The energies of the putative Dz^v global minima 
measured with respect to E avc , where E avc = — 2.605iV + 
2.5687V 2 / 3 - 4.859/V 1 / 3 + 8.614. 

is a ring of six interpenetrating icosahedra. The nearest- 
neighbour contact along the six-fold axis of this structure 
forms the common edge for six tetrahedra, whereas all 
other internal nearest-neighbour contacts are surrounded 
by five tetrahedra. This edge is said to have a discli- 
nation running along it,u and it is not possible to form 
a bulk compact polytetrahedral structure without in±iXfc. 
ducing such disclination lines. Frank-Kasper phasest§E£l 
and quasicrystals can, be characterized by their networks 
of disclination lines. □ Most of the clusters immediately 
above 7V=38 have structures based upon the Dz 38 global 
minimum.such as the 42- and 46-atom structures shown 
in Figure y. 

At larger sizes the global minima are rings of face- 
sharing icosahedra and double icosahedra that have a 
hole in the middle of the structure, e.g. 7V=52, 58, 60 
and 66 in Figure ||. Sometimes, as for N = 58, not all 
the icosahedral units are face-sharing, but may instead 
have a bridging tetrahedron between them. There are 
two forms for the 60-atom ring of six face-sharing icosa- 
hedra, which are analogous to the boat and chair forms 
of cyclohexane. The boat form has a slightly lower en- 
ergy (by 0.391) because a nearest-neighbour contact can 
be formed between the two icosahedra above the ring. 

At larger sizes three-dimensional network structures 
begin to form. For example, the 67-atom global minimum 
involves two perpendicular interlocking planar units, and 
the 86-atom structure is a three-dimensional network of 
eight, mainly face-sharing, icosahedra, with a complex 
pore structure (the six holes on the surface of the clus- 
ter connect to a single central pore). The formation of 
this latter three-dimensional structure leads to a substan- 
tial increase in stability (Figure ||). Above A^=86 similar 
structures dominate the lowest-energy minima that we 
have found. 

From these results we conclude that Dzugutov clusters 
tend to form non-compact polytetrahedral structures, 




FIG. 3: Some representative Dzjv global minima. The 
labels indicate the number of atoms in each cluster. 

based on packings of 13-atom icosahedra. Exceptions to 
this trend occur at N=b7 and 76, where the global min- 
ima are high-symmetry compact polytetrahedral struc- 
tures involving an ordered array of disclinations. As a 
consequence of this compactness these structures have 
a much higher strain energy than other global minima 
of similar size. They are also particularly stable for 
a modified Dzugutov potential that has been designed 
to favour compact-Bolytetrahedral clusters by having a 
lower value of n c g J?J and were only found on reoptimiza- 
tion of the structures obtained with this new potential. 
The structure oLihese new clusters will be described in 
detail elsewhere. L3 

Although the present work indicates that the lowest 
energy structures of small clusters are non-compact and 
polytetrahedral, a thorough study of the relative stabil- 
ities of possible bulk structures strongly suggests that 
for bulk-a body-centred cubic (bec) crystal is lowest in 
energy.tS Therefore, at sufficiently large sizes the lowest- 
energy clusters must be bcc. A simple comparison of 
-E avo , a polynomial fit to the energy of the global min- 
imum, with a similar polynomial function fitted to the 
energies of bcc rhombic dodecahedral clusters suggests 
that bcc clusters only become lower in energy beyond 
N~10 5 . 

For Dzugutov clusters, structures with disclination 
lines, although not dominant, can be competitive at some 
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TABLE II: Comparison of the DZ13 and LJ13 PESs. The samples of stationary points contain n m j n minima and 
nts transition states. For the n sea rch lowest-energy minima transition state searches were performed parallel and 
anti-parallel to the Hessian eigenvectors with the n cv lowest eigenvalues. Of the transition states found, were 
connected to the global minimum. (tWn) is the average geometric mean frequency of a minimum and (i^t™) is 
the average magnitude of the imaginary frequency of a transition state. AE, b u and ba are the energy difference, 
uphill barrier and downhill barrier, respectively, between two minima connected by a single transition state. S is the 
integrated length of the pathway between the two minima and D is the distance between—the two minima. N, the 
cooperativity index, gives a measure of the number of atoms involved in the rearrangement. olJ The average values of 
these quantities over all the non-degenerate rearrangement pathways are given below. (A degenerate pathway connects 
different permutational isomers of the same minimum.) (AE) = (b u ) — (bd)- The units of energy and frequency are e 
and (e/^eq) 1//2 i respectively, facilitating comparison with results for 13-atom Morse clusters.E3 
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sizes. Similarly for bulk, under certain conditions (when 
the system is somewhat compressed) a Frank-Kasper <j- 
phase is lowest in energy.Eil 

In recent simulations of a bulk Dzugutov supercooled 
liquid, domains corresponding to extended icosahedral 
configurations have been found to grow as the temperiu. 
ture is decreased and lower-energy minima are sampled. E£l 
Much lower diffusion rates are associated with atoms in 
these domains compared with atoms in the rest of the 
liquid. The configurations of these icosahedral domains 
look very similar to the types of non-compact polyte- 
trahedral clusters that we observe for clusters. Both the 
structures of the clusters and the domains have a common 
origin in the form of the potential. In particular, the po- 
tential maximum causes polytetrahedral structures to be 
favoured, and of the possible polytetrahedral structures 
the width of the potential favours non-compact forms 
based on icosahedra. 

Our results also indicate that high-density, compact 
polytetrahedral structures are disfavoured because they 
are associated with excessive values of -Entrain and E nnn . 
This result has implications for the form of the phase 
diagram, and in particular the stability of a dense liq- 
uid. Most simulations of Dzugutov liquids have been 
performed at constant volume. However, if a Dzugutov 
crystal is heated at a low constant pressure then there 
is a transition temperature at which the system jumps, 
straight to a low-density phase that resembles a foam.E2l 
Further simulations suggested that there is-*o region of 
the phase diagram where a liquid is stable£j 

Much work has been done on the effect of the range 
of a potential on the phase diagram. This research has 
shown that when the jC^p&c is sufficiently short there is 
no stable liquid phase.ocil The physical basis for this ef- 
fect is that as the range of the potential decreases, the 
potential well becomes narrower and so there is an in- 
creasing energetic penalty (i.e. -Entrain) for the dispersion 
of nearest-neighbour distances associated with the liquid. 
This destabilization of the liquid eventually leads to its 
disappearance from the phase diagramEJO 
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Roth attributed the absence of a liquid phase in t 
Dzugutov phase diagram to the range of the potential. 
However, the Morse potential with a similar width (i.e. 
p = p^ f f ) has a stable liquid phase. Only at much larger 
values of p (probably 13-14) does the liquid phase for the 
Morse potential become unstable. Instead, the present 
study clearly indicates that the absence of a dense liquid 
phase for the Dzugutov potential stems from the maxi- 
mum neart \/2r eq , because there are insufficient dense, 
non-crystalline structures with low values of both E nnn 

and -Entrain- 

Recently, Sear and Gelbart showed that adding a long- 
range repulsion of sufficient magnitude also leads to the. 
disappearance of the liquid from the phase diagram.efl 
As with the much shorter-ranged Dzugutov potential 
the positive value of E nnn destabilizes the liquid. Fu*-. 
thermore, as envisaged by Lebowitz and Pcnrose,E3 
Sear and GilbccU found that microphase separation in- 
stead occurred.E-3 Similar behaviour is seen for the 
Dzugutov potential: low-density non-compact polytetra- 
hedral structures are favourable for clusters, and Roth's 
heating simulations led to a low-density foam rather than 
a simple vapour.E3 Given this possibility of low-density 
disordered network phases, calculation of the true phase 
diagram for the Dzugutov potential is of obvious interest. 



B Energy Landscapes 

For simple pair potentials, such as the Morse and LJ 
forms, the topography of the -EES for 13-atom clusters 
has been examined in detail,Eili29 using near complete 
sets of minima and transition states. Here we examine 
some of the consequences of the Dzugutov potential for 
the topography of the energy landscape by comparing 
the Dzi 3 and LJ13 PESs. 

The properties of sets of stationary points for the two 
clusters are compared in Table 0. One of the main dif- 
ferences between the two systems is the much greater 
number of minima for DZ13. We obtained roughly twenty 
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times more minima for DZ13, yet the sample is far from 
complete. Transition state searches have been performed 
from less than 50% of the minima and far fewer searches 
were performed from each minimum than for LJ13. 

This greater number of minima is a consequence of the 
local maximum in the potential, which introduces barri- 
ers between configurations that would have been in the 
same basin of attraction for LJ13. In particular, there is 
now a barrier to the formation of nearest-neighbour con- 
tacts and so many Dz 13 minima involve atoms with low 
coordination numbers. For the LJ potential these atoms 
would be attracted by the long-range tail of the potential 
and so would increase their coordination. For example, 
when the DZ13 minima are reoptimized for the LJ po- 
tential 502 of them converge to the icosahedral global 
minimum. These differences are particularly evident in 
the connectivity of the low-energy regions of the PES. 
For LJ13 there are over 90 times more transition states 
connected to the global minimum. 

There are less dramatic differences between the other 
properties of the two PESs. The vibrational frequencies 
for the Dzugutov cluster are somewhat higher because 
of its slightly larger n c g value B3 For larger clusters this 
effect would be reinforced by the constraints on the sec- 
ond nearest-neighbour distances stemming from the max- 
imum in the potential. However, the latter effect is small 
for the 13-atom cluster because there are relatively few 
second nearest- neighbours. 

Barriers for the Dzugutov cluster are generally higher 
because the barrier for breaking a nearest-neighbour con- 
tact involves both loss of the contact energy and over- 
coming the potential maximum. In general, the barriers 
are also narrower and connect minima that are closer in 
configuration space, as expected for a PES with more 
minima. The rearrangements are also slightly less coop- 
erative. In summary, the DZ13 PES is much rougher than 
theLJ 13 PES. 

For larger clusters, there are generally a number of 
low-energy configurations associated with different ar- 
rangements of the icosahedral units that make up most 
low-energy clusters. Here we examine the effects of such 
competing structures on the PES topography, by calcu- 
lating the disconnectivity graph for a 37-atom Dzugutov 
cluster. For this cluster the global minimum, a linear 
quintuple icosahedron, is slightly lower in energy than a 
minimum based on the disk-like Dz3g global minimum 
with an atom removed. Next in energy above these two 
minima are a series of bent arrangements of five inter- 
penetrating icosahedra. 

Disconnectivity graphs provide ajpoeesentation of the 
barriers between minima on a PES E3c§Ell In a disconnec- 
tivity graph, each line ends at the energy of a minimum. 
At a series of equally-spaced energy levels we determine 
which (sets of) minima are connected by paths that never 
exceed that energy. Lines in the disconnectivity graph 
join at the energy level where the corresponding (sets 
of) minima first become corwteeted. In a disconnectivity 
graph an ideal single-funnec3'c3 PES is represented by 




FIG. 4: Disconnectivity graph for DZ37. Only branches 
leading to the 200 lowest-energy connected minima are 
included. The graph was calculated from a set of 30 000 
minima and 33 520 transition states. The five lowest- 
energy minima are depicted alongside the appropriate 
branches of the graph. 

a single dominant stem associated with the global min- 
imum to which the other minima join directly. For a 
multiple-funnel PES there are a number of major stems 
which only join at high energy. 

The DZ37 disconnectivity graph clearly has a multiple- 
funnel character. Each of the five lowest-energy struc- 
tures are separated by large barriers, and the correspond- 
ing branches all join up at the node at E=— 70.5. For 
example, the barriers for the lowest-energy pathway be- 
tween the funnels of the two lowest-energy minima are 
8.43e and 8.33e. These large barriers will strongly influ- 
ence the dynamics, making the rate of interfunnel pas- 
sage extremely slow. Hence, the low-energy minima act 
as efficient kinetic traps. 

The results for the two clusters in this section also 
show why the global minima of Dzugutov clusters axe. 
particularly difficult to locate. Bot. h . |PE S jpughnessEil 
and multiple-funnel PES topographyE3CJE3'E2l have been 
shown to hinder global optimization. The latter is partic- 
ularly troublesome, when, as with DZ37 , the global mini- 
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mum is at the bottom of a narrow funnel and is thermo- 
dynamically most stable only at low temperatures. 

These results also have implications for the dynamics of 
the supercooled Dzugutov liquid. Low-energy liquid cop, 
figurations with different extended icosahedral domainsEa 
will be separated by large energy barriers in a manner 
similar to DZ37. Therefore, as the average energy of the 
liquid minima decreases with temperature and the icosa- 
hedral domains grow, we expect the effective activation 
energy associated with diffusion or viscosity to inctease, 
thus leading to fragile (non-Arrhenius) dynamics S3cB 

C Conclusions 

The global minima of Dzugutov clusters have non- 
compact polytetrahedral structures. The interior atoms 
in these clusters almost always have an icosahedral co- 
ordination shell, rather than other polytetrahedral co- 
ordination shells that involve disclination lines, because 
the latter involve too large a strain energy. The 13-atom 
icosahedral units that make up the global minima are 
usually interpenetrating or face-sharing, and both these 
ways of linking icosahedra are commonly observed in 
metallic crystals that have icosahedral coordination.113 As 
the size of the clusters increases the arrangements of the 
icosahedral units change from one-dimensional chains to 
two-dimensional rings to three-dimensional porous net- 
works. These results are consistent with the previous ob- 
servation that the bulk liquid can sustain a distribution 
of free volume in the form of atomic size vacancies. E-3 

Unfortunately, Dzugutov clusters, because of their 
non-compact character, do not provide a realistic model 
for polytetrahedral metallic-clusters, such as those re- 
cently observed for cobalt. til However, promising re- 
sults have been obtained from a modified version of 
the Dzugutov potential that is better able to accom- 
modate the strain associated with compact polytetrahe- 
dral structures. Results for this system will be reported 
elsewhere. E3 

Dzugutov clusters typically have rough multiple- funnel 
energy landscapes. These features makes global opti- 
mization difficult and have implications for the dynamics 
of supercooled Dsugutov liquids, where extended icosa- 
hedral domains^ which are similar to the structures 
we find for the clusters, grow as the temperature is de- 
creased. Furthermore, our results help to explain why 
the bulk phase diagram appears not to include a stable 
liquid, and suggest the possibility of phases with complex 
low-density network structures. 
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